Computational engineering analysis of external geometrical modifications on MQ-1 unmanned combat aerial vehicle

This paper focuses on the effects of external geometrical modifications on the aerodynamic characteristics of the MQ-1 predator Unmanned Combat Aerial Vehicle (UCAV) using computational fluid dynamics. The investigations are performed for 16 flight conditions at an altitude of 7.6 km and at a constant speed of 56.32 m/s. Two models are analysed, namely the baseline model and the model with external geometrical modifications installed on it. Both the models are investigated for various angles of attack from −4° to 16°, angles of bank from 0° to 6° and angles of yaw from 0° to 4°. Due to the unavailability of any experimental (wind tunnel or flight test) data for this UCAV in the literature, a thorough verification of calculations process is presented to demonstrate confidence level in the numerical simulations. The analysis quantifies the loss of lift and increase in drag for the modified version of the MQ-1 predator UCAV along with the identification of stall conditions. Local improvement (in drag) of up to 96% has been obtained by relocating external modifications, whereas global drag force reduction of roughly 0.5% is observed. The effects of external geometrical modifications on the control surfaces indicate the blanking phenomenon and reduction in forces on the control surfaces that can reduce the aerodynamic performance of the UCAV.     

共轴刚性旋翼气动干扰数值计算方法

建立了一个适用于共轴刚性旋翼气动特性分析的数值模拟方法。该方法采用任意拉格朗日欧拉方法(Arbitrary Lagrange Euler,ALE)描述的可压缩Navier-Stokes(N-S)方程求解流场,采用低数值耗散的Roe格式进行空间离散;使用多重嵌套网格方法以模拟双旋翼的运动。针对共轴刚性旋翼配平,引入"差量修正"策略解决了传统配平中雅克比矩阵计算复杂的问题。首先,对Harrington-2共轴双旋翼的悬停气动性能进行了计算,然后,对某2 m直径共轴双旋翼的悬停及前飞状态进行了计算,并与试验值进行了对比。结果表明:在典型状态下拉力系数的计算结果与试验值误差在3%以内,扭矩系数的计算结果与试验值误差基本在5%以内;所采用的数值计算方法对旋翼涡尾迹特征具有较高的捕捉精度,可以有效模拟共轴刚性旋翼悬停和小速度前飞下的复杂流场及其细节特征。     

基于非结构/混合网格的高阶精度DG/FV混合方法研究进展

DG/FV混合方法因其具有紧致、易于推广获得高阶格式及相比同阶精度DG方法计算量、存储量小等优点,自提出以来已成功应用于一维、二维标量方程和Euler/N-S方程的求解。综述了DG/FV混合方法的研究进展,重点介绍了DG/FV混合方法的空间重构算法、针对RANS方程的求解方法、隐式时间离散格式、数值色散耗散及稳定性分析、计算量理论分析,并给出了系列粘性流算例的计算结果,包括用于验证混合方法数值精度的库埃特流,以及方腔流、亚声速剪切层、低速平板湍流、NACA0012翼型湍流绕流等。数值计算结果表明DG/FV混合方法达到了设计的精度阶,且相比同阶DG方法计算量减少约40%,而隐式方法能大幅提高定常流的收敛历程,较显式Runge-Kutta的收敛速度提高1~2个量级。     

旋转机翼悬停气动特性研究

鸭式旋转机翼(CRW)是一种先进的高速直升机方案,旋翼同时也是机翼,是其最关键气动部件之一。采用结构分区拼接网格技术进行空间离散,分区建立参考系,通过求解多重参考系下的N-S方程来计算旋翼流场,首先以传统的Caradonna-Tung实验旋翼的亚、跨声速悬停流场分析为例验证该方法的可靠性,进而采用该方法对旋转机翼悬停流场进行了数值计算,旋翼拉力计算值和地面实验值吻合较好,结果分析表明旋转机翼的悬停流场有着不同于传统旋翼的流场特性。     

战损飞机动力学特性分析及杀伤研究

现代战争中,防空武器性能不断提升,对飞机造成的威胁越来越大,确定飞机受到打击时的杀伤状态及杀伤等级成为飞机设计时需要考虑的重要问题。飞机被防空武器击中时,可能发生机翼脱离损伤,其质量、面积产生明显变化,将导致飞机的重心、气动特性的变化,严重危害飞机的稳定性。针对上述情况,以某型作战飞机为例,建立损伤情况下的飞机模型,综合分析机翼脱离损伤对飞机重量、重心和机翼面积带来的影响;建立损伤飞机的气动模型,计算损伤飞机气动特性的变化情况。基于这些变化,分析损伤对飞机动力学特性的影响,建立损伤状态下的动力学方程,计算损伤飞机的失控时间,依此确定飞机的杀伤等级。本文所做研究可用于在飞机设计过程中分析其杀伤状态及杀伤等级,为飞机设计时飞机杀伤的分析评估提供一种有效的解法。     

A coupling VWM/CFD/CSD method for rotor airload prediction

A coupling fluid-structure method with a combination of viscous wake model(VWM),computational fluid dynamics(CFD) and comprehensive structural dynamics(CSD) modules is developed in this paper for rotor unsteady airload prediction. The hybrid VWM/CFD solver is employed to model the nonlinear aerodynamic phenomena and complicated rotor wake dynamics;the moderate deflection beam theory is implemented to predict the blade structural deformation; the loose coupling strategy based on the ‘delt method' is used to couple the fluid and structure solvers.Several cases of Helishape 7A rotor are performed first to investigate the effect of elastic deformation on airloads. Then, two challenging forward flight conditions of UH-60 A helicopter rotor are investigated, and the simulated results of wake geometry, chordwise pressure distribution and sectional normal force show excellent agreement with available test data; a comparison with traditional CFD/CSD method is also presented to illustrate the efficiency of the developed method.     

三维可压缩 Navier-Stokes 方程的间断Galerkin 有限元方法研究

拓展了二维间断 Galerkin(DG)有限元方法研究，将该数值方法用于三维可压缩欧拉方程和 Navier-Stokes方程的求解。基于六面体网格单元，采用插值方法将物面的四边形面网格单元构造为弯曲面网格单元，更好地表述了真实物面特征；物面边界相邻体网格单元相应构造为高阶体网格单元，其余体网格单元采用八节点六面体单元，以较小的计算代价使网格满足 DG 方法计算需求。通过对三维带 bump 管道内流、圆球绕流以及旋转流线体绕流进行的数值求解，验证了边界弯曲方法的可行性及 DG 方法的高精度特性。此外，由于采用了隐式计算方法，仅需较少的时间步就能迭代收敛。     

网格拓扑对DLR-F6构型数值模拟的影响

基于TRIP3.0软件平台和多块结构网格技术,开展了控制方程和网格拓扑等因素对DLR-F6构型数值模拟的影响研究。数值计算采用与试验相同的参数,采用了RANS/TLNS方程,生成了O型和H型2种拓扑的粗、中、细网格进行模拟。分别从气动特性、压力系数分布曲线和表面流态3个方面对结果进行分析。通过与试验数据的对比表明,在模拟小分离流动时,采用H型拓扑网格和RANS方程模型,获得的计算结果会更为准确。     

Automatic differentiation based discrete adjoint method for aerodynamic design optimization on unstructured meshes

A systematic methodology for formulating, implementing, solving and verifying discrete adjoint of the compressible Reynolds-averaged Navier-Stokes (RANS) equations for aerodynamic design optimization on unstructured meshes is proposed. First, a general adjoint formulation is con-structed for the entire optimization problem, including parameterization, mesh deformation, flow solution and computation of the objective function, which is followed by detailed formulations of matrix-vector products arising in the adjoint model. According to this formulation, procedural components of implementing the required matrix-vector products are generated by means of auto-matic differentiation (AD) in a structured and modular manner. Furthermore, a duality-preserving iterative algorithm is employed to solve flow adjoint equations arising in the adjoint model, ensur-ing identical convergence rates for the tangent and the adjoint models. A three-step strategy is adopted to verify the adjoint computation. The proposed method has several remarkable features:the use of AD techniques avoids tedious and error-prone manual derivation and programming;duality is strictly preserved so that consistent and highly accurate discrete sensitivities can be obtained; and comparable efficiency to hand-coded implementation can be achieved. Upon the cur-rent discrete adjoint method, a gradient-based optimization framework has been developed and applied to a drag reduction problem.     

基于弹簧系统网格变形方法的旋翼气弹耦合分析

在旋翼气动弹性耦合（CFD/CSD）分析中引入弹簧系统网格变形方法，建立了一套适合于旋翼气动载荷分析的CFD/CSD耦合方法。为了解决CFD/CSD耦合中关键的网格变形问题，旋翼桨叶贴体网格变形采用基于“ball-vertex” 弹簧系统的动态网格方法，通过添加冗余约束，避免了畸形网格单元的产生。旋翼流场计算采用基于Navier-Stokes(N-S)方程的CFD模块，对基于运动嵌套网格 的空间流场进行求解，湍流模型采用B-L模型。结构分析采用基于中等变形梁理论的CSD模块，基于Hamilton变分原理建立旋翼桨叶动力学方程。首先对振荡NACA0012翼型的流场进行了求解，验证了网格变形模块和CFD模块的有效性，然后采用UH-60A直升机旋翼作为算例对结构动力学模块进行数值验证。在此基础上，计算了UH-60A直升机旋翼桨叶在前飞状态下的非定常气动载荷，并与飞行测试数据进行了对比。计算结果表明，文中的弹簧系统网格变形方 法可以有效地用于旋翼CFD/CSD耦合计算分析，提高了旋翼气弹载荷的预测精度。